A local fluctuation theorem for large systems 
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The fluctuation theorem characterizes the distribution of the dissipation in nonequilibrium systems 
and proves that the average dissipation will be positive. For a large system with no external source 
of fluctuation, fluctuations in properties will become unobservable and details of the fluctuation 
theorem are unable to be explored. In this letter, we consider such a situation and show how a 
fluctuation theorem can be obtained for a small open subsystem within the large system. We find 
that a correction term has to be added to the large system fluctuation theorem due to correlation of 
the subsystem with the surroundings. Its analytic expression can be derived provided some general 
assumptions are fulfilled, and its relevance it checked using numerical simulations. 
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INTRODUCTION 

In classical thermodynamics, nonequilibrium systems 
are irreversible. That is, the second law of thermody- 
namics stipulates that a macroscopic system evolves in 
one direction and fluxes have a definite sign. However, 
the second law is strictly applicable only to large systems 
or where properties are measured over long time scales. 
As the size of the system decreases, unusual events caused 
by thermal fluctuations become more frequent, and aver- 
age values of the flux with the opposite sign to that pre- 
dicted for a thermodynamic system are observable over 
finite periods. These effects are quantified by the Fluc- 
tuation Theorem (FT), which states that for a system 
initially at equilibrium and driven out of equilibrium, 
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where fit is an extensive function, the dissipation func- 
tion measured over a period of t |l| , which describes how 
irreversible a process is. In this equation, p(S7 t = A) 
refers to the probability that A takes on a value A ± dA. 
Eq. (JTJ is also asymptotically valid for steady-states. 
Initially based on a heuristic derivation and numerical 
evidence 0, fluctuation relations have been derived for 
closed systems in many different framewor ks (see [3, 0- 
5j for early results). Recently the fluctuation theorem 
has been used to measure physical quantities in experi- 
ments, for instance the torque of a molecular motor Q. 
In general, the use of the FT is relevant when the work 
done by the external forces is similar to the thermal en- 
ergy kgT (or its equivalent if another kind of thermal 
equilibrium is reached, e.g. 0). Otherwise, it reduces 
to stating irreversibility, because the probability of ob- 
serving a negative dissipation is so small that it can be 
considered impossible. Therefore, in large systems, (JTJ) 
can not be properly tested or applied. To deal with such 



a situation, one can record the dissipation over a small 
open subsystem where thermal fluctuations are sizeable. 
The dissipation in the subsystem will not satisfy equa- 
tion ([!} in general, and in this manuscript we consider if 
an analytical expression for a Local Fluctuation Theorem 
(LFT) can be obtained. 

Some experimental and theoretical work has been car- 
ried out on local fluctuation theorems since 1998. Mea- 
surements of local properties for steady states have been 
shown to satisfy fluctuation relations However, 
these studies employ an effective temperature, which can 
be considered to provide an ad hoc correction coefficient 
to the fluctuation relations. Gallavotti showed why the 
fluctuation theorem is expected to be valid at long times 
for the local entropy creation rate for weakly coupled 
systems Ayton et al. 113 considered the local fluc- 



tuation theorem for the dissipation function (JTJ , and pro- 
vided numerical results. In the present paper, we show 
that in a highly correlated system, the local dissipation 
obeys a LFT which is ([1]) with a linear correction term. 
The result explains the previous effective temperatures, 
the correction term can be analytically described in some 
cases, and it vanishes in the low correlation limit. 

In the next section we derive a LFT, then we numer- 
ically investigate the relevance of our assumptions in a 
realistic system. 



A LOCAL FLUCTUATION THEOREM 



Consider a large system that fulfils the assumptions 
of the FT, i.e. the invariance of the initial distribu- 
tion function under time-reversal mapping, ergodic con- 
sistency and time-reversibility of the equations of motion 
[l|. The FT can then be derived for an arbitrary phase 
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function <I>(r) odd under time-reversal \lc 
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If $ is chosen to be the dissipation function, @ becomes 
the usual FT. Here, we split S7 into two contributions Sle 
and where 0^ is the dissipation function measured 
over an arbitrary volume of length £, and choose $ = f^. 
In this case, becomes a local fluctuation theorem 
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This equation provides an exact expression for the cor- 
rection term. If there is no correlation between the two 
local dissipations, this term vanishes and £li obeys a bare 
FT (JTJ. Otherwise, some assumptions have to be made 
to obtain an effective description of this term. We will 
make two main assumptions. The first one is that one 
can consider the local dissipation as a random variable 
instead of a phase function so that, 
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where a is a coefficient that depends on the shape and 
size of the subvolume, and £ a random variable that is not 
correlated with fi^. Physically, the first term describes 
the very strong correlations that exist near the bound- 
aries between the two domains, and £ stands for the dis- 
sipation far away from the volume of interest, which does 
not depend on Qej- Quadratic and higher order terms 
are neglected in but do not seem relevant. This rela- 
tion between the two local dissipations leads to an exact 
expression for the LFT, 
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This model explains why the left hand side of (J5), called 
the asymmetry function, has been experimentally found 
to be a linear function of A with a slope that is not unity. 
Given our assumption of non-correlation between £ and 
Qi t , the coefficient a is 
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In solid states, typical lengths of correlation in the dissi- 
pation are typically much smaller than that of the sub- 
volume that has to be considered, and this correction 
term can be neglected. Significant exceptions are sys- 
tems close to a phase-transition, in which case an ana- 
lytic expression for a can often be derived (e.g. for the 
Ising model). On the other hand, in liquid and gas, cor- 
relation lengths are larger as particles can travel thought 
several subvolumes, and we shall see that they can be 
related to diffusion lengths. We focus on such systems, 



consisting a fluid driven out of equilibrium by an external 
field that produces a dissipative flux. This is very general 
and includes the studies on fluidized granular medium Q 
and Rayleigh-Benard convection but can also model 
Poiseuille or Couette flows, diffusion processes, and so on. 
To simplify the notation we consider the case where the 
field, F e and dissipative flux J are in the same direction 
and then the dissipation function reads 
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where j3 is the inverse temperature to which the system 
would relax in the absence of the field, and we adopt the 
notation that the average value of J is positive in a dis- 
sipative system. Then a refers to the spatial correlations 
of the integrated flux Jt : 
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Another expression for this coefficient can be derived af- 
ter defining details about the subvolume. We consider a 
rectangular unit cell of length C and with field applied in 
the x-direction. The subvolumes are obtained by divid- 
ing the cell into slices orthogonal to the direction of the 
field with width t (see figure [T]), but note that the follow- 
ing computations can be adapted to other situations. In 
this case, the relevant correlations arc fully described by 
the function C(x) defined by 
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where for Jg %t centred on x, jt(x) = lim^o Je,t/£ is the 
flux density at x. This function describes the decay of 
spatial correlations : if there is no correlation between 
the flux in these two sections, as when x goes to infinity, 
this term is equal to zero. On the other hand, C(0) = I 
due to the normalisation. A typical correlation length in 
the fluid can then be defined by 



C{x)dx 
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This length is to be compared to the subvolume 's to know 
whether a bare fluctuation theorem ([T]) is expected to be 
valid (£ £q) or if a LFT is required. 

We then make a second general assumption : noting 
that the function C satisfies C(0) = 1, C(oo) = and 
has a typical decay length of £q, we assume it is modelled 
by an exponential decay C(x) = e^ x ^°. Considering © 
this implies that (j t {0)j t {x)) - (jt(0))(jt{x)) = Ae~^°. 
a is then explicitly computable and assuming C is large, 
([5]) becomes 
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This LFT provides an analytic expression for the cor- 
rection term, which vanishes in the low correlation limit 
(£ ~> £q) and can be used instead of an effective temper- 
ature. In derivation of ([lip , we consider the large system 
limit, £-£^>£ Q . 

We now consider a particular case and show that it is 
possible to derive an expression for the correlation length 
under some conditions. Our system consists of N parti- 
cles of charge ci subject to a field in the x direction, and 
the dissipative flux becomes Jt = J Yli=i c i v x,ids where 
v Xj i is the x component of the velocity for particle i. For 
a system close to equilibrium, the correlation length, £q 
takes its origin in the Brownian motion of the particles, 
and for large C the following equalities hold : 
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where Ax = x(t) — x(0) is the x-displacement for one 
particle and C the length of the system, which does not 
appear in the final equilibrium expression for £q, 
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The last equality comes from the assumption of a Gaus- 
sian distribution for Ax. If the external field is not too 
high, £q is close to its equilibrium value and (|14|) provides 
a much more convenient way to measure it (directly or 
via the Einstein relation). 

Therefore, we have seen that for a fluid (or gas) driven 
out of equilibrium by an external field, the local dissipa- 
tion recorded in a section of length £ fulfils ([TT]) with £$ 
given by (fT4|) provided : the system is large, the field not 
too large (so that £q is close to its equilibrium value) and 
the decay of correlations is reasonably described by an 
exponential function. 



NUMERICAL STUDY 



where q i5 and c, = (—1)* are the coordinates, mo- 
menta and colour of the ith particle (cj = for wall 
particles) and Si is a switch equal to 1 for the wall parti- 
cles and zero otherwise. Fj is the interparticlc force on a 
particle, derived from a Weeks- Chandler- Anderson short- 
ranged repulsive pair potential Q , k the strength of the 
traps that fix the positions of wall particles, and F e the 
external field which induces a flux J = X)j=i c i v x,i- Fi- 
nally, a is a Gaussian thermostat that fixes the kinetic en- 
ergy of the walls : the fluid particles are not thermostated 
and obey their natural dynamics. The simulation used 
320 particles in a two-dimensional space with periodic 
boundary conditions. The wall temperature was set at 
1, the wall density at 0.8, the fluid's at 0.4, the field at 
F e = 0.08, the length of system at £ = 50.6 and the inte- 
gration time at t = 80. All trajectories started from the 
equilibrium distribution. The volume was divided into 
subvolumes as described above. This system is showed 
figffl 

The mean value of the total dissipation is (f2 t ) ~ 36 
and in 500 000 samples no negative value was observed. 
On the other hand, the local dissipation recorded in the 
volume of length £ showed in figfJl was 14.4 times smaller 
than the total volume. It consequently produced a num- 
ber of negative values of fit and its asymmetry function 
can be computed. As expected, it is a linear function 
with a slope larger than one, see fig. O 

The correlation length £q was measured out of equilib- 
rium by fitting Cix) = e~ x l l ° to the simulated values of 
C(x) and was found to be very close to its equilibrium 
value, related to a diffusion process and given by (fT4|) . cf. 
figlSl This tends to show that £q would often be experi- 
mentally equal to its equilibrium value, as external fields 
used in molecular dynamics simulations are very large 
(in this simulation, the mean velocity eventually reaches 
a few percent of the thermal one). The main assumption 
of the derivation was that an exponential decay e~ x ^° 
would fit the function C(x), which is not exact, cf. figlU 

Using the calculated value of £ = 4.56 at t = 80, 
and I — 3.52, the slope predicted by (fTTj) is 3.30 which 
is in good agreement with the numerically determined 



To present an application of this LFT and check the 
relevance of our assumptions, a system of color-charged 
particles between atomic walls was numerically studied. 
Color charges are similar to electric charges in that the 
particles feel a force proportional to their charge from 
the field, but there are no Coulomb interactions between 
the particles. The system is a high density gas in a long 
channel surrounded by thermostated walls. The dynam- 
ics are, 



q, = 



(15a) 




p 4 = Fj(q) + c t F e e x - S t {ap i + fc(q - q )) (15b) 
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FIG. 1. A schematic diagram of the system studied. Wall 
particles are thermostated and represented in black (a = 0), 
whereas fluid particles can be blue (a = 1) or red (a = —1). 
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FIG. 2. PDFs for the local dissipation function in a sub- 
volume 14.4 times smaller than the total volume and its Gaus- 
sian fit. Its mean value is {Qt,t) — 2.5 and the asymmetry 
function, pe(A) = In (p(Qe,t = A)/p(Qe.t = —A)), is found to 
be a straight line of slope 3.147 ± 0.001. 
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FIG. 3. Evolution of the correlation length £o determined out 
of equilibrium using a fit of C'(x) — e~ x ^° to the simulated 
data at various times. The dashed line is its theoretical equi- 
librium value determined using (|14p and the Einstein relation. 



slope of 3.146 ± 0.001 determined from the data shown 
in fig. [2] In figlSJ the predicted and actual slopes are 
compared for a range of /, and are shown to be in very 
good agreement. Therefore, even if the correlation decay 
is not exactly exponential, which could happen as this 
function is system-dependant, the predictions of this LFT 
remain robust. 
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FIG. 4. The decay of spatial correlations in the integrated flux 
calculated with ([9J (crosses) and is exponential model e~ x ^ l ° , 
where £o is the equilibrium correlation length estimated as- 
suming a diffusive process (solid line). 
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FIG. 5. Slope of the asymmetry function (left hand side of 
(|f f [) ) versus the length of the subsystem. The dashed line is 
the theoretical result derived in section 2 with the coefficient 
to estimated at equilibrium via a diffusive process. 



CONCLUSION 

In this letter, we considered a class of large systems 
(liquid or gas driven out of equilibrium by an external 
field). A local fluctuation theorem is derived for the dis- 
sipation measured in a section orthogonal to the field and 
provides an explicit expression for the slope of its asym- 
metry function. The proof can be adapted to describe 
other shapes of subvolumes. Therefore, arbitrary effec- 
tive temperatures no longer have to be used. This result 
holds provided some assumptions are fulfilled, mainly a 
low external field and an exponential decay of correla- 
tions in the flux. According to the numerical study, this 
last condition is not restrictive, and a reasonable accor- 
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dance is enough for this LFT to be verified. Finally, we 
pointed out that the need to resort to a LFT depends 
on the ratio £/£q, where I is the length of the subvol- 
ume over which the local dissipation is measured and £q 
is related to an equilibrium diffusion process and can be 
easily measured via the Einstein relation. As discussed 
in [13[, this FT is also asymptotically valid for steady 
states with an additional assumption of a form of decay 
of correlations in the dissipation with time. 
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